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Abstract 



VLBI is unique, among the space geodetic techniques, in its contribution to defining and main- 
taining the International Celestial Reference Frame, providing precise measurements of coordinates of 
extragalactic radiosources. The quest for increasing accuracy of VLBI geodetic products has lead to a 
deeper revision of all aspects that might introduce errors in the analysis. The departure of the observed 
sources from perfect, stable, compact and achromatic celestial targets falls within this category. This 
paper is concerned with the impact of unaccounted frequency-dependent position shifts of source cores 
■ in the analysis of dual-band S/X VLBI geodesy observations, and proposes a new method to measure 

' them. The multi-frequency phase transfer technique developed and demonstrated by Middelberg et al. 

, (2005) increases the high frequency coherence times of VLBI observations, using the observations at 

a lower frequency. Our proposed Source/frequency phase referencing method endows it with 
astrometric applications by adding a strategy to estimate the ionospheric contributions. Here we report 
on the first successful application to measure the core shift of the quasar 1038+528 A at S and X-bands, 
and validate the results by comparison with those from standard phase referencing techniques. In this 
particular case, and in general in the cm-wavelength regime, both methods are equivalent. Moreover 
the proposed method opens a new horizon with targets and fields suitable for high precision astrome- 
Qh trie studies with VLBI, especially at high frequencies where severe limitations imposed by the rapid 

^ , fluctuations in the troposphere prevent the use of standard phase referencing techniques. 

1. Introduction 

Geodetic VLBI observations with a network of antennas at the Earth are affected by the propaga- 
^ ■ tion medium, mainly the ionosphere and the troposphere. It is a basic practice in geodesy to cali- 
brate the ionospheric contribution with simultaneous observations at S/X-bands (2.2GHz/8.4GHz, 
respectively). The ionospheric- free delay observables at X-band (r^) are estimated from a combi- 
nation of observed delays (t x , t s ) at both bands: 



T'*-' X T"" & ST" 

and used to estimate the geodetic parameters. This approach works under the critical assumption 
that the brightness distributions for each source are identical and are co-located at both frequencies. 
In sources for which the VLBI core position is frequency dependent the exact expression must 
include a 24-hour sinusoidal extra term whose amplitude depends on the magnitude of the shift of 
the position of the source core (core shift): 



+ ^ * Ar|r ~ 0.08 * Arfr 
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where Arf^ = - sx , and 6 SX corresponds to the core shift between S and X-bands. 

The non-inclusion of this extra term introduces errors in the estimated ionosphere- free observables, 
and hence on the astrometric/geodetic products from the analysis. For sources with non- varying 
(i.e. stable) core shift 9 SX , the unaccounted extra term will propagate into an offset from the true 
X-band position, of magnitude (9 SX * 0.08) in the direction away from the S-band position [8]. 
Instead, unstable core shifts can propagate also into uncertainties in the estimated Earth orienta- 
tion parameters (Engelhardt, these proceedings) in the multi-epoch geodetic analysis. Section 2 is 
concerned with the origin, magnitude and geodetic impact of core shifts; Section 3 discusses the 
methods to measure them and Section 4 presents the results of our proposed method. 

2. Core shifts do exist 

Changes in the observed core positions at different frequencies have been measured in several 
sources, for example, 1038+528 A [4], [9], 4C39.25 [1], 3C395 [2] and 1823+568 [7]. We propose to 
classify the core shifts in two groups depending on their origin: 

• "Astronomical core shifts" result from opacity effects in the jet. The unresolved "core" of 
a compact extragalactic radio source is believed to mark the location where the optical depth to 
synchrotron self absorption ~ 1. This position changes with observing frequency as R CO re <x ^ _1//fcr , 
where k r depends on physical conditions in the jet. Core position shifts between S and X-bands 
of up to 1.5 milli-arcsecond (mas) are predicted in [3]. 

• "Instrumental core shifts" result from convolving the source structure with different resolutions 
at different frequencies, causing core shifts of up to half the beam size at the lower frequency. 

While the existence of source "core shifts" cannot be predicted, there are some clues which alert 
one to them. Larger "astronomical core shifts" are expected for flat spectrum sources, where the 
power index k r ~ 1; "instrumental core shifts" can be expected if the source structure at the 
higher frequency falls within a small fraction of the beam size at the lower frequency. Regardless 
of its nature, both core shifts have an identical effect on the analysis of S/X geodesy data. Table 
1 lists the propagation of plausible unaccounted stable core shifts into the analysis products. 



Observing 


"Astronomical" 


"Instrumental" 


Source position 


frequencies 


core shift 


core shift 


error 


2.2/8.4 GHz 


— 1 mas 


— 2 mas 


- 200 fias 


8.4/30 GHz 


- 0.25 mas 


— 0.5 mas 


— 50 /ias 



Table 1: Propagation of unaccounted non- varying core shifts into an offset from the true X and ri- 
band position in the geodetic analysis of 2.2/8.4 GHz (S/X) and 8.4/30 GHz (X/K) observations, 
respectively. 
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Figure 1. Hybrid map of 1038+528 A at S-band (left), and at X-band (right). The S-band beam super- 
imposed on the X-band map illustrates the case of structure blending effects from insufficient resolution at 
lower frequencies, and therefore "instrumental" core shifts. 

3. Ways to measure frequency-dependent core shifts 

Using closure phase relations in the image processing of VLBI data results in lack of absolute 
positional information in the hybrid maps. A rigorous alignment of maps at different frequencies, to 
measure frequency-dependent shifts of the core position, requires absolute astrometry observations, 
or standard phase-referencing to a nearby (achromatic) radio source. If astrometric observations 
are not feasible, a simpler but more imprecise procedure for extended sources is to use an optically 
thin component to align maps at different frequencies, and close epochs, and then estimate the 
change in the position of the core. 

Recently, Middelberg et al. [5, 6] proposed a new astrometric method that uses fast frequency 
switching observations of the target source and relies on the transfer of calibration from the lower 
to the higher frequency, after scaling by the frequency ratio. Their implementation proved to be 
a successful strategy to calibrate the rapid fluctuations of the troposphere, and hence extended 
the coherence time, in VLBA observations at 86 GHz, using interleaved scans at 15 GHz. This 
allowed the detection of a very weak, 100 mJy source. It also served to unveil the non-integer 
frequency ratio problem in the application of this method. On the other hand, the unaccounted 
dispersive ionospheric contamination, which was non-negligible even at these high frequencies, pre- 
vented them from making a proper astrometric measurement of the core shift. 

We present an extension of this method, a so-called "Source/frequency phase referencing" 
which complements the fast frequency switching observing strategy with source switching, in order 
to calibrate the remaining dispersive contaminating contributions to the observables. Section 4 
describes the first successful astrometric measurement of a core shift with this method, and Ap- 
pendix A contains a brief discussion of the basics of the method. 
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4. Source/frequency phase-referencing 

Conventional VLBI at high frequencies is severely constrained by the short coherence times 
imposed by the rapid fluctuations in the troposphere. The non-dispersive nature of the tropo- 
spheric propagation makes it possible to use lower frequency observations to calibrate higher, 
providing the switching interval between frequencies matches the temporal structure of the tro- 
pospheric fluctuations (i.e. coherence time) at the lower frequency. This is the basis of the 
multi- frequency phase transfer technique developed and demonstrated by Middelberg et al. [5, 6]. 
Our Source/frequency phase referencing method adds a source switching observing strat- 
egy, in addition to the fast frequency switching, to calibrate the dispersive contributions in the 
phase transfer strategy. The nodding between sources has to match the temporal and spatial 
structures of ionospheric propagation, and other non-dispersive terms, such as instrumental based 
contributions. A complete description of the method is given in Appendix A. 

We have successfully applied this method to the astrometric analysis of dual-band S/X VLBA 
observations of the pair of quasars 1038+528 A and B, 33" apart, and measured the core shift 
in the quasar 1038+528 A. The calibration transfer between frequencies involves multiplying the 
phases by the frequency ratio. The calibration transfer between sources is done as in standard 
phase referencing. The result is a Source/frequency phase referenced map, shown in figure 
2, whose offset from the center is a direct measure of the combined core shifts in the two sources 
between the two frequencies. The interpretation in terms of individual contributions from each 
source, assuming that shifts in the core position in each source would occur along the jet axis 
directions, is simplified by the quasi-orthogonal structures in this pair. The close alignment of 
the core shift offset with the A quasar source axis suggests a dominant contribution arising from 
this quasar; moreover, the quantitative agreement between the magnitude of the offset and the 
separation between the core and second component in the map of quasar A at X-band suggests 
the "instrumental" dominant nature of the offset. The results from a previous standard phase 
referencing analysis [9] are in complete agreement with those presented here, validating this new 
approach. 

5. Conclusions 

Unaccounted core shifts in the analysis of dual frequency VLBI geodesy observations propagate 
into offsets from the true X-band positions for the ICRF. We estimate deviations up to 200 /j,as 
and 50 /xas, respectively, for S/X and X/K observations, assuming temporally stable core shifts. 
Moreover, unstable core shifts can also corrupt the estimated Earth orientation parameters. 
We have successfully applied the Source/frequency phase referencing method to the ana- 
lysis of dual band S/X VLBA observations of the pair of quasars 1038+528 A and B, and measured 
a core shift in quasar 1038+528 A of ca. 800 fias. Our result is equivalent to, within the errors, 
to those obtained using standard phase referencing techniques [9]. As far as we know this is the 
first case of successful astrometric application of this multi- frequency phase transfer method. In 
this case the simultaneous observations of both frequencies and sources allowed solutions despite 
the non- integer ratio of the observed frequencies. While this new strategy does not present any 
advantage with respect to traditional techniques in the cm-wavelength regime, it does hold a big 
potential at high frequencies, which are out of the range of conventional phase-referencing. In par- 
ticular, we foresee a big impact when applied to observations of molecular line emission, where it 
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Figure 2. Left: Source/frequency phase referenced map of 1038+528 B from S/X VLBA observa- 
tions; the ~ 800/ias offset in NE direction is an estimate of the combined core shift of 1038+528 A and 
B quasars between S and X bands. Right: Same map, with the hybrid map of 1038+528 A quasar super- 
imposed, to show the agreement between the offset and the separation between the 2 components in the 
hybrid map of A quasar. This is an argument in favor of "instrumental core shift" . 

may provide bona fide astrometric alignment of emission arising from different transitions and help 
to elucidate the controversy between the proposed pumping mechanism for masers in evolved stars. 



Appendix A. The basics of the new method 

This section outlines the basics of this astrometric method Source/frequency phase referencing 
aimed to measure core shifts in radio sources. Its application involves observations of the target and 
a nearby source (in the formulae, A and B) at the two frequencies of interest (here x and s). At the 
post-processing, the amplitude calibration is done using traditional techniques, for all observations; a pure 
self-calibration analysis is used to solve for the phase, delay and rate of each low frequency (s) target source 
(A) observations. Following the standard nomenclature, the phase values <f> 8 A are shown as a compound of 
geometric, troposphcric, ionospheric and instrumental terms - assuming that structural contributions <\> a A str 
have been computed using the hybrid maps, and removed: 

<f>A = <t>A,geo + ^Itro + ^A.zon + <t>A,inst + 2 ™A > with n A inte g er 

These values are scaled by the frequency ratio, R, and used to calibrate the high frequency observations. 
The resultant high frequency referenced phases to the low frequency are: 

+ ^Xinst R^hnst) + Mn A ~ R.n s A ) (1) 

This calibration strategy results in perfect cancellation of non-dispersive troposphcric terms: 
i^Xtro - R 4Xtro) = . but not for thc dispersive ionosphere: {4>\ ion - R-<f>A, ion ) = (R~ ^Wa^u 
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Taking this into account in Equation (1) above results in: 

ft R4 S A = <P\str + 2ttz,^4^ + (R- i)^ iion + (d>* A . nst ~ R4 s AMst ) + Mn A R-n s A ) (2) 

where 2-kv D ' 6 *' sx — (<j>A, ge o ~ R4 S A,geo)' an( i ®A,sx is the core shift in A between s and x. 

Similarly, the analysis of the observations of a nearby calibrator, B, after removing structural terms, <f) B atr 
and 4> S B str at both frequencies, results in: 

</>% - R4% = 2^^%^ + (R- ±W B , lon + (<f> B ,in S t - R4 S B,inst) + M^B ~ R.n' B ) (3) 

which are transfered for further calibration of A source observations. The resultant Source/frequency 
referenced phases, combining (2) and (3), are: 

(ft - R-ft) - {4>% -R4%) = ft.str + 2 ^- ( ^T fe ' 3x) +(R-Tt) (<P S A,ion ~ <^,o„) + i<f>A. ms t ~ R-ft.inJ 

- i^BMst ~ R4 s B , mst ) + M{n A - R.n' A ) - K - R.n%)] 

A careful planning of the observations, namely switching between sufficiently nearby sources with a duty 
cycle which matches the ionospheric/instrumental time-scale variations, would result in negligible differen- 
tial ionospheric error and instrumental corruption, that is: 

(#-s)(^-^ o J~0 

{^Xinst ~ R4 S A,inst) ~ {<t>B,inst ~ R 4h,inst) ~ 

which results in an expression for the Source/frequency referenced phases for the target source free 
of ionospheric/instrumental corruption: 

{ft-R4 A )-((j) B -R4 B ) = (p x A str +2-KV D - {eA ' sx ~ eB ' sx) +2-Kn" , withn" integer if R integer (ornj, = n%) 

And finally, the calibrated complex visibilities from the target observations are inverted to yield a synthesis 
image of A at x-band, where the offset from the center is an estimate of the combined core shifts in A and 
B between s and cc-bands. 
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